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ABSTRACT 

Millisecond radio pulsars are neutron stars that have been spun-up by the transfer of angular momen- 
tum during the low-mass X-ray binary phase. The transition from an accretion-powered to a rotation- 
powered pulsar takes place on evolutionary timescales at the end of the accretion process, however it 
may also occur sporadically in systems undergoing transient X-ray activity. We have obtained the first 
optical spectrum of the low mass transient X-ray pulsar SAX J1808. 4-3658 in quiescence. Similar to 
the black widow millisecond pulsar B1957-I-20, this X-ray pulsar shows a large optical modulation at 
the orbital period due to an irradiated companion star. Using the brightness of the companion star as 
a bolometer, we conclude that a very high irradiating luminosity, a factor of ~ 100 larger than directly 
observed, must be present in the system. This most likely derives from a rotation-powered neutron star 
that resumes activity during quiescence. 

Subject headings: accretion, accretion disks — binaries: close — stars: neutron 



1. INTRODUCTION 

SAX J1808. 4-3658 was the first-discovered low mass X- 
ray binary (LMXB) transient showing coherent pulsations. 
This confirmed unambiguously the long-sought connec- 
tion between LMXBs and millisecond radio pulsars (Bhat- 
tacharya & van den Hevuel 1991; Tauris & van den Hevuel 
2004). The detection of coherent X-ray pulsations during 
outbursts testifies that the neutron star possesses a mag- 
netic field oi B ^ 10^ — 10^ G, sufficient for a small mag- 
netosphere to form (Psaltis & Chakrabarty 1999; Menna 
et al. 2003; Di Salvo & Burderi 2003). Unlike persistent 
LMXBs, SAX J1808.4-3658 is a transient system, i.e. it 
is active in X-rays only for short intervals lasting a few 
months (outbursts) followed by quiescent periods of years. 

During quiescence LMXB transients are very faint in 
X-rays (5-6 orders of magnitude less than in outburst) 
usually with luminosities of 10^^ — 10'^^ erg s^^. Tran- 
sient systems, therefore, represent a unique laboratory for 
the study of compact objects in accretion regimes that are 
inaccessible to persistent sources. Most neutron star tran- 
sients are characterized by a quiescent X-ray spectrum 
consisting of a soft component, usually ascribed to the 
cooling neutron star which has been heated during out- 
bursts (Brown et al. 1998; Rutledge et al. 1999), and 
a hard power-law tail (with photon index F ~ 1 — 2) 
the nature of which is still debated (Campana & Stella 
2000). Several mechanisms have been put forward to ex- 
plain these spectral components, ranging from accretion 
disks (in different flavours such as advection-dominated or 
convection-dominated disks, disks stopping at the magne- 
tosphere, etc.; Narayan et al. 1997; Blandford & Begel- 
man 1999; Igumenshchev et al. 2003) to emission from 



the interaction between the relativistic wind from a re- 
activated radio pulsar with matter outflowing from the 
companion star (Stella et al. 1994; Campana et al. 1998b). 
XMM-Newton carried out the first detailed study of SAX 
J1808.4-3658 in quiescence: SAX J1808.4-3658 was dim- 
mer than any other neutron star transient (with a 0.5-10 
keV luminosity of 5 x lO'^^ erg s~^ at a distance of 2.5 kpc, 
in't Zand et al. 2001) and exhibited only a hard power- 
law component (F = 1.4^q 3, 90% confidence level) in its 
quiescent X-ray spectrum (Campana et al. 2002). These 
results were confirmed by a subsequent XMM-Newton ob- 
servation (Wijnands 2002). 

In the optical, SAX J1808. 4-3658 is dim during quies- 
cence (mean R ~ 20.9±0.1; Homer et al. 2001) and bright- 
ens considerably in outbursts (i? ~ 16.2 ± 0.2; Wang et al. 
2001). The optical light curve in outburst and quiescence 
is modulated at the orbital period and it is in anti-phase 
with the X-ray light curve, likely indicating that irradi- 
ation of the companion star plays a crucial role in spite 
of the low X-ray luminosity. This is unlike other quies- 
cent transients. The mass function derived from X-ray 
data (4 x 10"^ Mg, Chakrabarty & Morgan 1998) and the 
requirement that the companion fills its Roche lobe led 
to the conclusion that it must be a rather low mass star, 
possibly a brown dwarf (Bildsten & Chakabarty 2001). A 
white dwarf companion is ruled out because it would not 
fill its Roche lobe. Homer et al. (2001) proposed that the 
bulk of the optical emission in quiescence arises from the 
internal energy release of a remnant disc and the orbital 
modulation from the varying contribution of the heated 
face of the companion star. Burderi et al. (2003) noted 
that the required irradiating luminosity needed to match 
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the optical flux, however, is a factor 10 — 50 higher than the 
quiescent X-ray luminosity of SAX J1808.4-3658. These 
authors proposed an alternative scenario, in which the irra- 
diation is due to the rotational energy emitted in the form 
of a relativistic particle wind from the fast spinning neu- 
tron star, which switched to the rotation-powered regime 
during quiescence. Their results are in agreement with the 
weak constraints from the optical magnitudes by Homer 
et al. (2001). 

2. DATA ANALYSIS 

Here we present indirect evidence for an active radio pul- 
sar in SAX J1808. 4-3658 during quiescence with the first 
optical spectra and / band photometry of SAX J1808.4- 
3658. These data were obtained with the ESO-VLT (UT4 
Yepun) during two half nights on July 12-13 2002. We 
carried out / band photometry with 3 min exposures with 
F0RS2 (pixel size of 0".126/pixel and a field of view of 
6'. 8 X 6'. 8) over one orbital period. Spectroscopy of the 
same target was performed using the low resolution grism 
600RI (centered at 6780 A with a resolution of 55 A/mm) 
covering 5120 8450 A and a 1" slit with 3 min spectra over 
four orbital periods. Data reduction was done in MIDAS 
to remove the bias level, flat-field. 

The region around the optical companion of SAX 
J1808. 4-3658 is crowded and poor seeing conditions com- 
plicated the analysis (varying between 1.5 — 2"). We take 
advantage of previous ESO-VLT images (obtained during 
quiescence in 1999 with seeing ~ 0.5") to de-blend our data 
(these magnitudes were V = 21.82±0.03, R = 21.63±0.04 
and / = 21.08 ± 0.04). Our I photometry shows a dimmer 
source (/ = 21.5 ± 0.1) and a clear modulation at the 2.01 
hr orbital period with a semi-amplitude 0.2 ± 0.04 mag 
(65% in flux. Fig. 2). We also re-calculated the modula- 
tion semi-amplitude in V (0.13 ±0.06) and R (0.39 ±0.09) 
using the Homer ct al. (2001) data. The /-band light 
curve shows a clear maximum at phase 0.52 ± 0.05 and a 
single minimum at phase 0.02 ± 0.05 (based on the precise 
X-ray ephemerides Chakrabarty & Morgan 1998). This 
is a clear indication of emission from an irradiated com- 
panion (e.g. Charles & Coe 2004) and it argues against 
emission from the impact point between the gas stream 
from the companion and an accretion disk (the hot spot) 
since this has maximum at phase 0.8 — 0.9. 

We obtained spectra at 3-min intervals over four or- 
bital periods. We selected spectra taken with seeing better 
than 1.6" (due to poor seeing conditions), collecting a to- 
tal of 51 min of good data. Wavelength calibrations used 
HeArNe arc lamp observations. Second-order flexure ef- 
fects were corrected using night-sky emission lines. This 
correction was always < 0.3 A. Spectra were corrected for 
slit-losses according to Diego (1985). We also account for 
the contaminating stars, estimating their relative contri- 
bution in a 1" slit on the good-seeing VLT images and 
interpolating to the spectral range. Errors were tracked 
along these processes resulting in a 0.1 mag error. A 
weak Ha emission line is visible in the spectrum (equiva- 
lent width EW = 10.3 ± 3.7 A, 68% confidence level, and 
FWHM = 44.0 ± 6.3 A, see insert in Fig. 1). 

3. MODELLING THE DATA 



We first checked that at the time of our observations 
the source was in quiescence. A Rossi X-ray Timing Ex- 
plorer (RXTE) pointed observation, performed 15 days 
before the optical observations, provided a 3it upper limit 
of ^ 10'^^ erg s"-'^ (0.5-10 keV). Similar limits are provided 
by the Galactic bulge scan carried out with RXTE (near- 
est observation on Jul 10, 2002; Swank et al. 2002). WhUe 
these limits are not particularly constraining, they are suf- 
ficient to exclude that the source was in outburst (note 
that an outburst started three months later Wijnands et 
al. 2003). 

Since SAX J1808. 4-3658 was in quiescence, what is the 
cause of the optical emission? Likely candidates are emis- 
sion from the companion star and/or the disk. In order 
to fit within the Roche lobe of a 2.01 hr binary the com- 
panion mass has to be less than 0.17 M©. In the model of 
Bildsten & Chakrabarty (2001) the most likely companion 
is a 0.05 Mq brown dwarf bloated by irradiation to fill 
its Roche lobe (0.13 Rq). The maximum intrinsic opti- 
cal luminosity from the companion for any of the models 
by Bildsten & Chakrabarty (2001) is ~ 3 x lO^^ erg s'^ 
(corresponding to a star temperature of 4800 K for a dis- 
tance of 2.5 kpc). This is too low a luminosity to account 
for the observed optical flux, which is a factor of > 10 
brighter. We therefore turn to the accretion disk as a 
possible source. Assuming that the quiescent X-ray lumi- 
nosity is powered by accretion we can infer the expected 
mass inflow rate and derive the; corresponding optical lu- 
minosity (including irradiation), which fails to account for 
what we see by more than a factor of 100. Disk models 
may be envisaged with a much higher mass accretion rate 
together with a truncation radius (fine) tuned to avoid op- 
tical and soft X-ray violation of observed data. However, 
this kind of models still require some additional ingredient 
to explain the large optical phase modulation. The emis- 
sion from the piflsar could itself extend to the optical and 
extrapolation of the power law X-ray fiux (assuming the 
XMM-Newton observation found the X-rays in a similar 
state) could account for about 30% of the optical luminos- 
ity. But this also could not explain the observed orbital 
phase modulation. 

It is instructive to compare the properties of SAX 
J1808. 4-3658 with those of the black widow pulsar PSR 
B1957+20 (Fruchter et al. 1988), which consists of a 
1.6 ms radio pulsar irradiating its white dwarf compan- 
ion (orbital period of 9.16 hr) with a rotational energy of 
]^q35 gj-g g-i g^j-^j 7-rays are generated in an inter- 
binary shock front, which causes ablation and heating of 
the companion (Phynney et al. 1988; Arons & Tavani 
1994). An X-ray nebula has recently been revealed around 
PSR B1957±20 confirming this scenario (Stappers et al. 
2003), and the orbital modulation is large with an i?-band 
semi-amplitude > 4 mag (Callanan et al. 1995). A sim- 
ilar system is the eclipsing millisecond radio pulsar PSR 
J2051-0827, consisting of a 4.5 ms pulsar orbiting its very 
low mass companion (~ 0.03 Mq) every 2.4 hr (Stappers 
et al. 1996). Radio eclipses as well as a ^ 3.3 mag op- 
tical modulation have been observed in PSR J2051-0827 
(Stappers et al. 2001), however no X-ray observations are 
available. 

Inspired by this analogy and following Burderi et al. 
(2003), we now attempt to account for the optical and 
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X-ray spectra (even if not close in time) as well as the 
V, R and / modulations of SAX J1808.4-3658 with an 
irradiated star plus the contribution of the shock front. 
We fit the data by using the irradiating luminosity (Ljrr), 
the fractional luminosity difference between the heated 
and the cold face of the companion (/) and interstellar 
absorption {Ay) as free parameters (e.g. Chakrabarty 
1998). We obtain a good fit to all the available data 
(reduced = 0.7 with 57 degrees of freedom'', see Fig. 
1). In particular, the required irradiating luminosity is 
Li„ = [4+1) X 10^3 erg s'^ (90% confidence level for 
three free parameters, i.e. Ax^ = 6.3). The best fit 
fraction is / = 0.65 ± 0.10 resulting in temperature dif- 
ference at the two faces of the companion star of about 
1000 ±300 K. This temperature difference is similar to the 
one observed in PSR J2051-0827. The amplitude mod- 
ulation of SAX J1808. 4-3658 is instead smaller than in 
the case of PSR J2051-0827 and PSR B1957-h20. This 
could be due to some remaining contamination due to our 
de-blending process, or it could be an inclination effect. 
For the most strongly modulated PSR J1957+20, the in- 
clination is 70°. PSR J2051-0827 has an inchnation of 
~ 40°. 

The estimated absorption is Ay = 1.0 ±0.5, in line with 
previous estimates. Optical emission from the shock front 
accounts for about 15% of the total emission**. In the fit we 
assumed that all the irradiating luminosity is re-emitted 
by the star. Therefore the irradiating luminosity we de- 
rived represents only a lower limit since a relativistic par- 
ticle wind could have rather different effects from those of 
X ray irradiation in terms of effective albedo. Moreover, 
we note that evidence of this large irradiating luminos- 
ity comes from the equivalent width of the Ha line. If 
the line comes from reprocessing, one can roughly infer 
an irradiating luminosity of ^ 3 x 10^'^ erg s~^. This 
can be estimated by taking the line flux and increasing 



it by the fraction of emitted flux intercepted by the star 
{R^/Aa?, with i?, the companion star radius taken equal 
to the Roche lobe radius and a the orbital separation) and 
by the fraction of energy re-emitted in Ha (e.g. Hynes et 
al. 2002). Here we assumed a conservative value of 0.3 
(see the discussion in Hynes et al. 2002). 

4. CONCLUSIONS 

The required irradiating luminosity is in all cases large 
(~ 4 X 10^3 erg s~^), indeed much larger than the ob- 
served X-ray luminosity in quiescence; neither accretion- 
driven X rays, nor the intrinsic luminosity of the com- 
panion star or disk are able to account for it (see also 
Burderi et al. 2003). The only source of energy avail- 
able within the system is then the rotational energy of 
the neutron star. In order to have such a large spin- 
down luminosity one needs a neutron star magnetic field 
5 ^ 6 X lO'' G (the companion star albedo might be larger 
than zero and the neutron star emission may also be par- 
tially beamed; this value is smaller than the one estimated 
by Burderi et al. 2003, since they did not account for the 
contamination from nearby stars). The required magnetic 
field is well in the range inferred from X-ray observations 
during the outbursts (see above). Prospects for directly 
observing SAX J1808.4 3658 pulsing in the X ray band 
are rather hard since with XMM- Newton we collected less 
than 300 counts in 30 ks (Campana et al. 2002). Given the 
crowding around SAX J1808. 4-3658 in the optical band, it 
would be difficult to detect an Ha nebula like that around 
PSR B1957+20. On the other hand a search in the ra- 
dio band looking for a millisecond radio pulsar would re- 
quire searches at high frequencies to overcome the effects 
of free-free absorption (Campana et al. 1998a; Ergma & 
Antipova 1999; Burgay et al. 2003; Burderi et al. 2003) 
and a favorable orientation of the radio beam. 
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The conversion efficiency of rotational energy into X— rays for SAX J 1808.4—3658 is ^ 10 ^, a factor of 10 higher than PSR 1957+20 (Stappers 
et al. 2003), as expected due to geometrical reasons. 

^ We fit our data with an irradiated star plus the extrapolation to the optical of the X-ray tail. The presence of a disk is not required by our 
fit. We tried in any case to fit with an irradiated star plus a disk model finding similar results for the star parameters and a large disk inner 
radius ( ^ 10® cm). 
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Fig. 1. Central panel: X~ray to optical (VLT) spectrum of SAX J1808. 4-3658. The data are corrected for absorption 
(Ay = 1.0) using standard extinction curve (Fitzpatrick 1999). Dots shows the XMM-Newton (X-rays) and VLT (optical 
spectrum plus / photometry) data. The continuous (blue) line shows the overall best fit for an irradiating luminosity 
of 4 X 10'^'^ erg s~^. The (green) dotted line shows the contribution from shock emission between the relativistic pulsar 
wind and matter outflowing from the companion. The power law tail has a photon index of F = 2.0, this provides a 
slightly better fit to the data with respect to a F = 1.5 model (in that case the required irradiating luminosity is higher). 
The lower continuous (red) line shows the mean contribution from the bloated (0.13 Rq) companion star. This spectrum 
together with the optical V , R and / modulation has been fitted with an irradiated star plus shock emission model. Free 
parameters were the ratio between the 'hot' and the 'cold' star surface and absorption. A downhill simplex method for 
the search of the minimum has been adopted. Errors were computed in the 3-dimcnsion space with A^^ = 6.3. 
Upper insert: Optical spectrum with the best fit model spectrum (upper continuous [blue] line). The mean star spectrum 
(lower continuous [red] line) together with the minimum and maximum spectra (dotted [green] lines) are also shown. 
Lower insert: particular of the broad {FWHM = 44.0 ±6.3 A) Ha hue. The line equivalent width is EW = 10.3 ± 3.7 A. 
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Fig. 2. VLT / band light curve (lower panel) plus V (upper panel) and R (middle panel) taken from Homer et al. (2001) 
deconvolved light curves of SAX J1808. 4-3658. Filled dots represent the data and the blue-cyan heavy line is the best-fit 
model for the V data, the orange line the model for the R data and the red line for the / data. Phases are from the X-ray 
ephemeris, i.e. the mean orbital longitude as derived from the X-ray pulse arrival time delays. 



